Couplings between holographic dark energy and dark matter 
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We consider the interaction between dark matter and dark energy in the framework 
of holographic dark energy, and propose a natural and physically plausible form of 
interaction, in which the interacting term is proportional to the product of the powers 
of the dark matter and dark energy densities. We investigate the cosmic evolution 
in such models. The impact of the coupling on the dark matter and dark energy 
components may be asymmetric. While the dark energy decouples from the dark 
matter at late time, just as other components of the cosmic fluid become decoupled 
as the universe expands, interestingly, the dark matter may actually become coupled 
to the dark energy at late time. We shall call such a phenomenon as incoupling. 
We use the latest type la supernovae data from the SCP team, baryon acoustics 
oscillation data from SDSS and 2dF surveys, and the position of the first peak of the 
CMB angular power spectrum to constrain the model. We find that the interaction 
term which is proportional to the first power product of the dark energy and dark 
matter densities gives an excellent fit to the current data. 



I. INTRODUCTION 

A number of astronomical observations, such as type la supernovae (SNIa) [1], the cosmic 
microwave background (CMB) anisotropy and the large scale structure (LSS) [3|], have 
shown that the expansion of our universe is accelerating. The acceleration of the Universe 
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strongly indicates the existence of a mysterious exotic matter, named dark energy, which 
has negative pressure and takes the largest proportion of the total density in the Universe. 
The ultimate fate of our Universe is to be determined by the dark energy, yet we know 
little about its properties. Furthermore, while both the dark matter and the dark energy 
is commonly thought to interact weakly with ordinary matter, it is not immediately clear 
whether there is some interaction between the two, or if the two are in some way linked. 

In this paper we consider the dark matter-dark energy interaction in the framework of 
the holographic dark energy (HDE) model #tZ|- The basic ideas of the HDE model are 
based on the holographic principle [8[ , which was inspired by the Bekenstein entropy bound 
of black holes 2fl. In this model, the dark energy is related the quantum zero-point energy 
density p\, the total energy of the whole system with size L should not exceed the mass of 
a black hole of the same size, thus we have L 3 p^ < LM 2 V The largest infrared cut-off L is 
chosen by saturating the inequality so that the HDE density is 

p de = 3c 2 M*L- 2 , (1) 



where c is a numerical constant, and M p \ = l/vo7rG is the reduced Planck mass. If we 



take L as the size of the current Universe, for instance the Hubble scale H 1 , then the dar 
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subsequently 



energy density estimated with Eq.flTJ will be close to the observed value. £ 
pointed out that this yields a wrong equation of state for dark energy. Li 
proposed that the future event horizon might be chosen as the IR cut-off L, and viable HDE 
models can then be constructed. As demonstrated in Ref. the HDE scenario may help 
to answer the fine tuning problem and the coincidence problem of dark energy. 

Various aspects of the HDE model and its extensions have been investigated subsequently. 
For example, it has been generalized to the case of non-flat Universe [10|], and it was studied 
from the perspective of the anthropic principle (ljj. The holographic dark energy was 
also considered for driving the inflation, and its resulting power spectrum index calculated 
1121 1 . Ref. 13 1 defined a new future event horizon as the infrared cut-off L and found 
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that the coincidence problem could be alleviated to some extent. Ref. [14|] investigated 



the possible form of the cut-off scale L for the solution of coincidence problem. Inspired 



by the HDE model, Ref. [15[ proposed the Ricci scalar curvature radius as L. The model 
was reconstructed with a scalar field [if]]. Observational constraints have also been studied 
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The growth of density perturbations in the holographic dark energy models has 
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been studied in 



19]. 



Interactions between the DE and DM components were considered for the HDE model 
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2 If . Ref. 20] suggested that for the L = H 1 case, with interactions between the 



two dark components, the dark energy could acquire an equation of state in agreement with 
observations, thus partly solving the no-acceleration problem noted in Ref. [7] for interaction- 



free models. Ref. 



21] generalized the discussion in Ref. 20] to the case of a future event 



horizon. In that paper the proposed form of interaction is 3Hp tot , which could also lead to 
a negative dark energy equation of state. 

However, Ref. [22] pointed out that the interaction terms such as Q ~ 3Hp totl 3Hpd e , 
3Hp m all lead to unphysical results, because the dark matter density could be driven to a 
negative value in these models. 

Here we propose a form of interaction which we considered natural and physically plau- 
sible. In general, when two types of matter interact with each other, the interaction should 
be proportional to the product of the powers of the number density of both. This is most 
easily understood in the following way: if one component does not exist, then there cer- 
tainly would be no interaction at all; with its amount increasing, the interaction rate should 
increase. Thus, at least in the diffuse limit, the interaction rate should be proportional to 
some powers of the density (the relation could become more complicated when the densi- 
ties become high). For instance, this is the case for ordinary two-body chemical reactions, 
namely 

a" 3 ^ ~ n in2 („> , (2) 

i.e. the reaction rate is proportional to the product of number densities of the two particle 
species. 

However, besides the interaction in the usual sense ("scattering"), one may also consider 
a more general case where the dark energy is converted to dark matter or vice versa by 
"annihilation" or "decay". In such a case, the "interaction term" or more precisely the 
energy exchange term would depend only on one of the densities and not on the other. 

In general, we consider the interaction term to be proportional to some powers of the 
density, 

Pm + 3Hp m = "/PmPAi (3) 
p A + 3H(l + w A )p A = - 1P lpl (4) 
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in which p m and p\ are the energy densities of the matter and HDE respectively, and 7 is 
a parameter with mass dimension [density a+ ^ _1 xtime] _1 . The power index a, ft is assumed 
to satisfy 

a,p>0. 

where if one of them is 0, we will have the case of "decay" (if the other is 1) or "annihilation" 
(if the other is 2). If 7 > 0, the interaction suggests that dark energy is converted into dark 
matter, while 7 < suggests the inverse process. This form of interaction was considered 
in the study of coupled quintessence model 23J]. However, not much research on this form 
of interaction was carried out to understand how such a coupling term would affect the 
evolution of the Universe. 

In the next section, we present the evolution equations for the interacting holographic 
dark energy model (details of our derivation are given in the appendix). In Sec. 3, we study 
the dynamical behavior of the model. In Sec. 4 we explore the parameter space and find the 
best-fit and allowed region of a, (3. Finally, concluding remarks are given in Sec. 5. 



II. THE MODEL 



Now let us consider a spatially flat Friedmann- Robert son- Walker (FRW) Universe filled 
with matter component p m , HDE component p\ and the radiation p r . For simplicity we do 
not distinguish dark matter and ordinary matter below but incorporate both in the matter 
component p m . It is certainly possible or even likely that dark energy only couples to dark 
matter. However, in that case the baryons simply evolve as pb ~ a -3 where a is the scale 
factor, the baryon component would then become dynamically irrelevant as the Universe 
expands, so it only affects the details of the model but not its qualitative behavior, hence 
we will neglect this effect in the present paper. The Friedmann equation reads 

3M^H 2 = p m + PA + p T , (5) 

where p r is the energy density of the radiation which satisfies p r ~ a -4 . Thus, Eqs. (Jl]) , (J3j) , 
(H]), ()5]) are the dynamic equations which uniquely determine the evolution of p m (z), p\(z), 
Pr{z)i w\(z) and H(z). In Eq.((T]), L could be Hubble horizon H" 1 , particle horizon and 
future event horizon, i.e. 

L = R ph (t) = a(t) [ (6) 
Jo ti a z 
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or 



(?) 



L = R eh (t) = a(t) j 

J a 

We define the following dimensionless quantities: p m = p m /po cr , p r = Pr/Pocr, Pa = Pa/Pocv, 
where po cr = SHqM^ is the current critical density, and E(z) = H(z)/Hq, and introduce 



da 
~W7P' 

- Pm/ POcr, Pr 



A = 3 a+ ^Ho 



2(a+/3)-3 M 2(a+P)-2 
1V1 pl 



a+B-1 Tj-l 



(8) 



The dimensionless A parameter determines the characteristic strength of the interaction. 
The equations of motion of the HDE model for the particle (event) horizon cases are 



df>A 
dz 

dE(z) 
dz 

w A (z) 



Pa 



+ 
1 
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(9) 



in which the upper(lower) sign represents the case of particle (future event) horizon (c.f. the 
appendix for details of derivation). It is easy to verify that if L in Eq.([T]) equals the particle 
horizon R p h or Hubble horizon H , the model cannot fit the current observational data, 
because current data prefer the dark energy with equation of state near —1, while in the 
above two cases the equation of state w will be always greater than — | (see also the analysis 
in Ref. j3|). Therefore, in the following we will focus on the interacting model with L = R e h 
(see Eq.©). 



III. DYNAMICAL BEHAVIORS 



In this section we study the dynamical behavior in models with interaction between DE 
and DM, which may give rise to some interesting behaviors. First we consider the case where 
a, (3 are fixed to be 1, then consider the more general cases. 

As a typical example we consider the model with c = 0.81, A = —0.03, which happens 
to be the a = (3 = 1 model which best-fits the current data ( for details of the fitting 
procedure see next section). First we plot the evolution of densities of different components 
in FigfTJ While one might think that the coupling between DE and DM might help resolve 
the coincidence problem, we find that for these particular parameter values the evolution of 
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FIG. 1: Evolution of densities in the HDE model. The black-dashed line is the constant magnitude 
of dark energy in the LCDM universe. The evolutions of matter and radiation in LCDM model 
are very close to those (green-dashed and blue-dotted-dashed lines) in the HDE model. 

different densities is rather similar to that of the standard ACDM model, i.e. the DM and 
DE densities could still differ by orders of magnitude, so it does not seem to help much in 
this case. Rather, as in other HDE models, the fine tuning problem could be alleviated: the 
density of the dark energy is associated with the scale of horizon rather than the Planck 
scale, due to the holographic principle. We have also investigated cases with larger coupling 
constant A. We find that indeed, for the strong coupling case (A > 1), the dark matter 
and dark energy could be made to evolve with comparable density for an extended range of 
redshifts, thus in some sense solving the "coincidence problem". This is perhaps also true 
for any strongly coupled dark matter-dark energy model. However, for strong coupling the 
fit to observation is not good. 

To investigate how the dynamics of the components are affected by coupling with the 
other component, one may compare the expansion term 2>Hpi (i = m, A) with the interaction 
term 7p m pA- If the interaction term is greater than or comparable with the expansion term 
of both of the components, the two components are strongly coupled and evolve together. 
As the Universe expands, typically the density of the component would drop as a~ 3 ( 1+w \ 
the interaction term often drops much more rapidly than the expansion term. When the 
expansion term of a component is much greater than the interaction term, that component 
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FIG. 2: Comparisons of the expansion term (red solid line) for dark energy (left panel) and dark 
matter (right) with the interaction term(blue dashed line) for the model with c = 0.81, A = —0.03. 

evolves almost freely, and we may say the two components are decoupled 1 . The coupling 
and decoupling effects between the two components could be asymmetric, as each has its 
own expansion term. It is possible that for one component the expansion term is greater, so 
it is decoupled from the other component, but for the other component, its expansion term 
is still smaller than or comparable with the interaction, thus it is still coupled to the other 
component. 

In Fig. El we compare the expansion term of the dark energy and the dark matter with 
the interaction term as a function of redshift. In the left panel, we plot the magnitude of 
the dark energy expansion term 3(1 + w)p\, and in the right panel, we plot the magnitude 
of the dark matter expansion term 3p m , as red solid lines. The interaction term | A | p ™{^ 
is plotted in both panels for comparison with the expansion terms. In the early Universe, 
the dark matter takes a greater proportion of the total density, and for dark energy the 
interaction term is comparable with the expansion term, so the dark energy is moderately 
coupled to the dark matter. Later, as the density of dark matter is diluted by expansion, 
the interaction term becomes negligible compared with the expansion term, and the dark 



1 Note here that in the literatures the word "decoupling" may refer to a number of different phenomena. 
For the decoupling of WIMP dark matter from the radiation-baryon fluid, what is discussed above (the 
evolution of density) is usually called "thermal decoupling" . There is also "kinetic decoupling" , which 



refers to decoupling of temperature (see, e.g., 24|). Admittedly, this nomenclature is somewhat confusing, 
since one might naively think that "thermal decoupling" refers to decoupling of temperature, and "kinetic 
decoupling" refers to decoupling of density evolution! Below we shall refer to the decouplings discussed 
here as "dynamical decoupling", the terms "thermal decoupling" and "kinetic decoupling" are not used 
in this paper. 
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FIG. 3: Comparisons of the expansion term (red solid line) for dark energy (left panel) and dark 
matter (right) with the interaction term(blue dashed line) for the model with c = 0.2, A = 5. 

energy is decoupled from dark matter. In the right panel of Fig. [2j we plot the dark matter 
expansion term 3p m and interaction term A%^f as a function of redshift z. We see that 
throughout the cosmic history, for the dark matter the expansion term is always greater 
than the interaction term, so the dark matter is at most weakly coupled to dark energy. 
However, at lower redshifts, the difference becomes larger, and the dark matter become 
further decoupled from dark energy. Such a decoupling history is quite typical for many 
interactions. 

However, there could also be some quite unexpected and interesting behavior in the 
presence of dark energy. In Fig. [31 we make plots similar to Fig. [2] but for the model with 
c = 0.2, A = 5, and this model is chosen to illustrate such an unusual behavior. In this model, 
because w < — 1, the dark energy density actually increases at later time. We see from the 
left panel of Fig. [3j since initially the dark matter density is greater, that the interaction 
term is comparable to the expansion term for dark energy. However, later on the dark energy 
becomes dominant, and the interaction term falls below the expansion term of dark energy, 
so we may say that dark energy becomes decoupled. However, more interestingly, from the 
right panel of the figure, we see that initially the interaction term of the dark matter is 
smaller than the expansion term, so the dark matter is only weakly coupled to dark energy, 
but at later time, the interaction term is comparable with the expansion term of dark matter, 
so the DM becomes coupled to DE. This is the reverse of the familiar decoupling process, 
which we shall call incoupling. This strange and interesting behavior could only occur when 
a component is coupled to dark energy, whose density does not decrease during the cosmic 



9 



expansion. As the DE keeps being converted into DM, the DM evolution is greatly affected 
by DE. This model does not yield a good fit to the observational data, but is only used to 
illustrate an interesting dynamical behavior in the presence of w < — 1 dark energy. 

In FigJU we plot the evolution of H for three choices of c parameter, i.e. the best fit 
Co = 0.81, and cq ± cr c , while keeping A = —0.03. The Hubble parameter for the ACDM 
model is also plotted for comparison. We also plot the equation of state of the dark energy 
for these models in Fig. |5j On the left panel, curves for different c but identical A are plotted; 
on the right panel, curves for identical c but different A are plotted. In all of these models, 
the phantom divide w = — 1 is either crossed, or will be crossed in the near future. The 
dark energy density would keep increasing, the interaction with DM would not slow it down, 
as the dark energy has become decoupled to the DM. Eventually a "Big Rip" would occur. 
The current value of the equation of state Wq for the best-fit model c = 0.81, A = —0.03 
is around —1.04, which is within the la CL region of Wq in the WMAP 5-year results 251 ]. 
suggesting that this model is consistent with the current observations. 

The value of A in these models does not significantly affect the fate of the Universe, but it 
does greatly affect the past evolution. This result is easy to understand; as the dark energy 
dynamically decouples from dark matter, the interaction become negligible, regardless the 
precise value of coupling constant. The eventual fate of the Universe in this interacting 
holographic dark energy model is qualitatively similar to those without interactions, though 
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FIG. 4: Evolution of Hubble rate for three choices of parameter c, three curves (red, blue, and 
green one) are the best-fit, and ±1<t c C.L. values. Here we set f2 m o = 0.26,-f/o = 72km/s/Mpc. 
The black-dashed line is the Hubble parameter in the LCDM universe. 
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FIG. 5: Equation of state with a = (5 = 1 for several parameters. The black-dashed line is the 
constant equation of state in the LCDM universe. 

the details of the evolution would be different. Thus, as in the case of non-interacting 
holographic models, for the fate of the Universe, c < 1 leads to phantom-like behavior and 



Big Rip 26], while c > 1 leads to quintessence behavior. This is very different from the 
models discussed in Ref. For further discussions on the fate of the Universe and Big 
Rip singularity in our model, see Appendix B. 

We now consider the more general case of a, (3 ^ 1. We plot the equation of state in Fig. 
IH1 varying a in the left panel while keeping (3 = 1, and varying (3 in the right panel while 
keeping a = 1. The other parameters are kept as before. For a < 1, the equation of state for 
dark energy is slightly more negative. For a > 1, at higher redshift p^ increases rapidly, and 
the coupling term would become very large, the whole system would be dominated by the 
evolution of matter. However, such a model would not really fit the low redshift observation, 
so it would not be relevant to our purpose. For large (3, we can also see that the equation 
of state of dark energy w > — 1 at high z, but the variation is not as rapid as in the case 
of a, and this is because the dark energy component does not increase as much as the dark 
matter component. Still, this makes a poor fit to the observational data. 



Anyway, for the case a, {3 1, Eq.( 1A16l) is reduced to 



Mx) * -l§^(E{xf -Pa- p r ) a . (10) 

Substituting this equation into Eq. ()A4j) and making use of Eq. dAip . we obtain the following 
approximate equation for p~A 

Pa - 3p A , (11) 
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FIG. 6: Redshift evolution of the equation of state, with the best-fit values of c, A and f2 m o- The 
black-dashed line is the constant equation of state in the LCDM universe. 

for which the solution will simply be 

p A (a) ~ B/a\ (12) 
where B is a constant. This scaling behavior o 



dark energy density (if a, /3 1) is very 
similar to that of the Chaplygin gas model 27| ; which at the high redshift also has the 
scaling behavior as a -3 while at the low redshift has the constant magnitude. However, in 
the present case it does not provide a good fit to the observations. 



IV. FITTING THE MODEL 



Now we fit the model parameters with the current observational data. We use three 
kinds of observations; the first one is type la supernovae (SNe la) which serves as the 
standard candle; the second is the baryon acoustic oscillation (BAO) in large scale structure 
of galaxy distribution, the third is the cosmic microwave background anisotropics, the latter 
two provide standard rulers at different redshifts. For the flat geometry we are considering, 
d,L = (1 + z) 2 (Ia = (1 + z)r, where di,dA,r are the luminosity distance, angular diameter 
distance and co moving coordinate distance, respectively, and r is given by r = f* -^tj, 
which could be calculated for each model by solving Eq. ((9]) or Eqs. ( 1A17l) -( lAl~8l) . 

For the SNe la, we use the data set published by Supernova Cosmology Project (SCP) 
team recently 28j]. This data set contains 307 selected SNe la which includes several current 
widely used SNe la data sets, such as those collected by the Hubble Space Telescope (HST) 
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FIG. 7: Contours of A versus c for a = (3 = 1. 
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30|, the SuperNova Legacy Survey (SNLS) |3l| and the Equation of State: Sup Er Novae 



trace Cosmic Expansion (ESSENCE) surveys j^l]. Using the same analysis procedure and 
improved selection approach, all of the sub data sets are analyzed to obtain a consistent and 
high-quality "Union' SNe la data set, which provides tighter and more reliable constraints 
for our models. For the BAO data, we use the data given in Percival et al. {33 1 and WMAP 



team 



25]. We use the quantity r s /D v which is determined by the Sloan Digital Sky Survey 



(SDSS) (z = 0.35) and the Two Degree Field Galaxy Redshift Survey (2dFGRS) (z = 0.2) 
in our constraints, where r s is the comoving sound horizon size at the decoupling epoch, 



and D v is the effective distance defined by Eisenstein et al. [34] . Since this quantity is 
just a ratio of two "standard rulers" , it could give a reliable measurement of the expansion 
history of the Universe. For the CMB, we use the position of the first peak of the Cosmic 
Microwave Background (CMB) angular pow er spectrum t\. The details of this technique 



were described in our earlier paper Ref. 



35|. 



We use the Markov Chain Monte Carlo (MCMC) method to constrain the parameters 
with the combined data set of SN Ia+BAO+CMB (for details of our implementation of 
this technique, see Ref. jsoj] ) . In the present model, the relevant parameters are the matter 
density f2 m0 , Hubble constant H , and the holographic dark energy parameter c, A. We 
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FIG. 8: Probability distribution function of the power index parameters a and /3. 

marginalize over Hq and Q m o- For a,/3, we limit their range to < a, (3 < 2. As discussed 
in the last section, greater values of a, /3 would not yield a good fit to observation. 

For the fixed a = (3 = 1 case, as can be seen in Fig. [7J the best-fit value and error on c, A 

are 

c = 0.8ll{j;g, A = -0.03 ±0.02. (13) 

The behavior of the best fit model has been discussed as the example given in the previous 
section. The models within the la contour also exhibit similar behavior. 

Now we let a, (3 be free parameters, the PDF of a and are shown in Fig. [HI The best-fit 
values for a, (3 are a = 0.69^38 and (3 = 0.98^o; 8 ?. If one demand a, (3 to be integers for 
physical reasons, then a = 1, /3 = 1 is still the best fit. Alternatively, a non- integer power 
may rise if the interaction is complicated. The best fit for c and A are c = 0.8lion an d 
A = — 0.04+q;q4, as shown in Fig. [91 it is not too different from the a = (3 = 1 case. 

V. CONCLUSION 

The holographic dark energy (HDE) model is a dynamical model of dark energy, moti- 
vated by the holographic principle, in which the event horizon of the space-time plays an 
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FIG. 9: Contours of A versus c for free a and (5. 

important role in the dynamics. The microphysical basis of the model is not well under- 
stood, but the dynamics of cosmic expansion could be studied. Observationally, the HDE 
model is distinguishable from the ACDM, as its equation of state varies with time. The 
model could also be distinguished from other dynamical dark energy models. For example, 
in the HDE model the dark energy behaves as the Chaplygin gas at high redshift if a, /3 3> 1 
(see the discussion at the end of Sec. 3), but for the low redshift region, the equation of state 
of the Chaplypin gas approaches a constant, while for this model it continues to evolve. In 
properly constructed HDE models such as the ones discussed in this paper, the dark energy 
component is subdominant at high redshifts compared with matter or radiation, so it would 
not significantly affect the CMB or large scale structure power spectrum, though at low 
redshift it could affect the growth of structure in detail. We plan to investigate the details 
of the structure growth in subsequent studies. 

In the framework of the holographic dark energy model, we proposed a form of interaction 
between dark energy and dark matter which we considered to be natural and physically 
plausible. That is, the interaction rate is proportional to the product of powers of the dark 
energy and dark matter densities p^Pa- For this type of interaction, we find that for the weak 
coupling case (the dimensionless coupling constant A < 1, see Eq. (jHJ)), the ultimate fate of 
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the Universe is similar to the case of non-interacting holographic dark energy, namely for 
c < 1 the equation of state of the dark energy will cross w — — 1, and behave like a phantom 
in the future, and the Universe ends in a Big Rip. We also discovered the interesting 
phenomenon of incoupling, which occurs only when the dark matter is coupled to dark 
energy This is the opposite of the familiar process of decoupling, the dark matter becomes 
increasingly coupled to the dark energy as Universe expands. We then use the currently 
available observational data to constrain the parameters of this model. The data set we 
used include 307 selected SNe la published by the SCP team [281 ]. baryon acoustic oscillation 
from SDSS and 2dFGRS 33], as well as the position of the first peak of the CMB angular 
power spectrum 25[]. We find that the best-fit form of interaction is a = 0.69, = 0.98. If 
one requires a, (3 to be integers for reasons of fundamental physics, then the most favored 
interaction rate is proportional to the product of the density of dark matter and dark energy. 
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Appendix A: dynamical equations 

We give the detailed derivation of the dynamical equations of the interacting holographic 
dark energy model here. Using the dimensionless notations given in Sec. 2, we have 

E(xf = p m + p A + p r , (Al) 

p r = fi r0 e 4x , (A2) 

L - 3p m = -A§4' (A3) 

p;-3(l + WA )p A = A§4> (A4) 
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pa = ^ = 4r^ x)2 - (as) 

So 

3 4 

£(x) • E{x)' = -\p m + (1 + w A )p A + -p r ]. (A7) 

By substituting Eq. flA6P into Eq. fl A4[) . we have 

WA(X) - 1 3L(x) 3E(x)- (A8) 
If L is taken to be the particle horizon (Eq.Q) or future event horizon (Eq.((7])), the deriva- 
tive of L is 

L' = -(L±1), (A9) 

where the upper (lower) sign represents particle horizon (future event horizon). By sub- 
stituting Eq. ( 1A9|) and ( 1A5|) into Eq. ( 1A8j) . we have the equation of state in the particle 
horizon (upper sign) and future event horizon (lower sign) cases 

^ = ^l4n~l Mx) If ri < A1 °> 

3 3c-E(x) 3 E(x) 
Substituting Eq. ( 1A10I) into (1A7I) . we have 

3 

2E(x) ■ E\x) = 3(p m + ^p r ) + 2p A ± - ^rrP^A- ( A11 ) 

Taking the derivative of Eq. ( IA1I) with respect to x, we have 

2E(x)-E'(x)=p m + p A + 4p r} (A12) 
Substituting Eq. (1A3|) into Eq. (1A12[) . we obtain 

2£(x) • = 3p m + p A + 4p r - -?—~p a m p{. (A13) 



Subtracting equations ( 1A13I) from ( 1A13I) . we obtain the following simple dynamical equation: 



^w^i- < A14 > 



Eliminate p m in Eq. (1A11I) using Eq.t lAlj) . 



?, 



*v = § B <*> + 2k) <* ~ M ~ 2^ {E(xf ~ h ~ M " A ± ^krr (A15) 
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and 



= 4 ± sdfcy - \w) (Eix? - h - w * (A16) 



To make Eq. (1A14j) and (1A15[) easy to solve numerically, we change the variable to 



z. 



3 



dE(z) 1 r 3 „. , 1 . pi 



dz (l + z) L 2 w 2£(^) vr ' rJV c-£ 2 (^) 

A ;(E(z) 2 -p r -p A rff A ], 



2E(z) 2 

E(Q) = 1, (A18) 



where 



p r = a r0 (l + zf = 4.78 x 1(T 5 (1 + z) 4 , (A19) 

and the upper (lower) sign in Eq. (IA17I) and (IA18I) corresponds to the model with particle 
(future event) horizon. Thus, Eq. (1A17|) and (IA180 are the complete system of equations. 
Eqs. (IA16l) - flA18p reproduce Eqs. (jSJ). In the limit A = and p r = 0, the equations reduce to 
the non-interacting holographic dark energy model [3, 0]- The evolution of event horizon 
size L and the corresponding dark energy density are shown in Fig. [101 As z — > —1, the 
size of the event horizon shrinks to zero, and the dark energy density diverges, i.e. we are 
approaching a Big Rip singularity. 



Appendix B: Big Rip Singularity 



There have been extensive discussions on the future singularity of the dark energy models, 



particular 



Ref. 



38. 



y for the case where the future Universe is dominated by the phantom dark energy. 
391 ] investigated the appearance of the future singularity, and showed that the 



?he Big Rip is however, the only relevant 
39, W. Ref 



411 ] showed that the coupling 



singularities could be classified into four types, 
singularity in the holographic dark energy model 
between dark matter and dark energy may help to solve the coincidence problem, but not 
the Big Rip (Type I) singularity issue. In particular, [40] discussed the singularity problem 
in the holographic dark energy model, and proved that: 
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1. if holographic dark energy c = 1, there exists a de Sitter solution to the dynamic 
equation system of the HDE. But if c ^ 1, the de Sitter solution does not exist. 

2. For the de Sitter solution (c = 1), one should consider the back reaction of the quantum 
effects near the singularity. The contribution of the conformal anomaly as a back 
reaction near singularity could indeed moderate or even prevent the future singularity 
(Eqs.(53)-(61) in [40]). 

In our interacting HDE model, the best-fit parameter c = 0.81, and for such model the 
fate of the Universe is a Big Rip in the future. In the following, we investigate whether this 
conclusion is affected by the back reaction of quantum effects. 

Since the future of the Universe is dominated by dark energy, we will neglect all of the 
terms which contain p m and p r , as they will become subdominant when compared with the 
dark energy density p A . The dynamical equations then simplify to 

E(x) 2 =p A , (Bl) 

p' A -3(l + w A )p A = (B2) 

L' = -(L-±), (B3) 

where we only consider the case for future event horizon. From Eqs. ( 1B1I) -( 1B3|) . the asymp- 
totical behavior of H, p A and w A can be obtained easily: 

«A— (B4) 

PA = PA a 2( ^ 1} , (B5) 

H = H n% o a-°-\ (B6) 

where p A0 = SHqM^Q^q, and Q A0 is the fractional HDE density taken at the present day. 
From Eqs. flB6j) . the derivative terms are 

H= (H ni ) 2 (--Aa 2 ^, (B7) 



c 



i \ 3 fl 



2 



H {H nl j 2i:-l) a 3 ^\ (BS) 
H= (tf fi|o) 4 6Q-l) a 4( ^ 1} . (B9) 
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The conformal anomaly T A (back reaction term) has the following form 40| : 
T A = -12bH 2 + 2Ab'(-H 2 + H 2 H + H A ) - (46 + 66") (fi + 7HH + AH 2 + 12H 2 H). (BIO) 
Thus in the HDE model 



T A = l2Hfo 2 m f(b,c)a^- l \ 



(Bll) 



where 



/(ft, c) = 2b' 



-b 



1 


-36" 












[: 


!'- 




e-oi 


c 













1 



i-ll +7I--1I+4 



(B12) 



With Eqs. (1Blll) and ( 1B12I) . the energy density and pressure of the conformal anomaly p A 
are [40] . 

,4(1-1) 



i r* 

Pa = -- dta 4 HT A ~ -3Hfo 2 A0 cf(b,c)a 4 

a Jto 

PA = \ (T A + Pa) = 4 (l - |) H& 2 A0 f(b, c)a^\ 



3 v ' ' V I 
Under this back reaction, the expansion rate is 



(B13) 
(B14) 



H 



3^ 



(Pa + Pa) 



H 2 a 2 ^ 



n 



AO 



pi 



ttlof( b , c ) a 



2(^-1) 



(B15) 



Since the energy density of back reaction is positive (p A > 0), the function f(b,c) < 0; 
therefore, if 2 increase monotonically in the future if c < 1. Thus, we prove that for the c < 1 
case, the back reaction effect from quantum corrections cannot prohibit the occurrence of 
the Big Rip. 
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